Precise measurement of the positive muon anomalous magnetic moment 
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A precise measurement of the anomalous g value, = [g — 2) /2, for the positive muon has been 
made at the Brookhaven Alternating Gradient Synchrotron. The result = 11 659 202(14)(6) x 
10"^" (1.3 ppm) is in good agreement with previous measurements and has an error one third 
that of the combined previous data. The current theoretical value from the standard model is 
a;,(SM)= 11 659 159.6(6.7) x 10'^° (0.57 ppm) and a^,(exp)-a^(SM) = 43(16) x 10"^° in which 
a;j(exp) is the world average experimental value. 
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Precise measurement of the anomalous g value, 
= {g—2)/2, of the muon provides a sensitive test of 
the standard model of particle physics and new informa- 
tion on speculative theories beyond it. Compared to the 
electron, the muon g value is more sensitive to standard 
model extensions, typically by a factor of (m^/me)^. In 
this Letter we report a measurement of for the posi- 
tive muon from Brookhaven AGS experiment 821, based 
on data collected in 1999. 

The principle of the experiment, previous results, and 
many experimental details have been given in earlier pub- 
lications 0J^. Briefly, highly polarized /Li"*" of 3.09 GeV/c 
from a secondary beamline are injected through a su- 
perconducting inflector ||^ into a storage ring 14.2 m in 
diameter with an effective circular aperture 9 cm in diam- 
eter. The superferric storage ring Q has a homogeneous 
magnetic field of 1.45 T, which is measured by an NMR 
system relative to the free proton NMR frequency 
Electrostatic quadrupoles provide vertical focusing. A 
pulsed magnetic kicker gives a 10 mrad deflection which 
places the muons onto stored orbits. The muons start in 
50 ns bunches and debunch with a decay time of about 
20 /zs due to their 0.6% momentum spread. Positrons are 
detected using 24 lead/scintillating fiber electromagnetic 
calorimeters read out by waveform digitizers. The 
waveform digitizer and NMR clocks were phase-locked 
to the Loran C frequency signal. 

The muon spin precesses faster than its momentum 



rotates by an angular frequency uja in the magnetic field 
{B) weighted over the muon distribution in space and 
time. The quantity is 



where oja is unaffected by the electrostatic field for muons 
with 7 = 29.3. Parity violation in the decay /x"*" — > 
e^iJ^Vf, causes positrons to be emitted with an angular 
and energy asymmetry. Because of the Lorentz boost, 
the positron emission angle with respect to the muon 
spin direction in the muon rest frame is strongly coupled 
to its energy in the laboratory frame. The number of 
decay positrons with energy greater than E is described 
by 

N{t) = 7Vo(S)e-*/(''^) [1 + A{E) sm{uat + ME))] (2) 

in which the time dilated lifetime 7r w 64.4 /is. Some 
140 g — 2 periods of 4.37 /is were observed. 

Most experimental aspects of the data taking in 1999 
were the same as in 1998 [|l|. However, some improve- 
ments were made. Care was taken in tuning the AGS 
ejection system to minimize background from any extra- 
neous proton beam extracted during the muon storage 
time. Scintillating fiber detectors which could be moved 
in and out of the storage region were used to study beam 
properties. Scintillation counters in front of five calorime- 
ters were used to measure muon losses from the storage 
ring. 
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The principal new feature of the 1999 run is the 20- 
fold increase in data collected. During a data acquisition 
time of 500 hours, we obtained about 2.9 billion decay 
positrons with energies greater than 1.0 GeV. The AGS 
delivered 4 x 10^'^ protons of 24 GeV in six bunches, sep- 
arated by 33 ms, over its 2.5 s cycle. This resulted in 
5 X 10'* stored muons per cycle. 

The magnetic field B is obtained from NMR measure- 
ments of the free proton resonance frequency ujp. Sev- 
enteen NMR probes are mounted in an array on a trol- 
ley which moves on a fixed track inside the muon stor- 
age ring vacuum. The trolley probes are calibrated with 
respect to a standard spherical H2O probe to an accu- 
racy of 0.2 ppm before and after data-taking periods. 
Interpolation of the field in the periods between trolley 
measurements, which are made on average every three 
days, is based on the readings of about 150 fixed NMR 
probes distributed around the ring in the top and bottom 
walls of the vacuum chamber. Figure |l| shows a mag- 
netic field profile averaged over azimuth. The variations 
in the amplitudes of the multipoles affect (B) by less 
than 0.02 ppm. The average readings of 36 uniformly 
distributed fixed probes are maintained to 0.1 ppm by 
feedback to the main magnet power supply. 
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FIG. 1. A 2-dimensional multipole expansion of the field 
averaged over azimuth from one out of seventeen trolley mea- 
surements. One ppm contours are shown with respect to a 
central azimuthal average field Bo = 1.451 266 T. The cir- 
cle indicates the muon beam storage region. The multipole 
amplitudes at the beam aperture radius of 4.5 cm are given. 

The 1.45 T magnetic field of the superconducting in- 
fiector is well shielded from the storage region by flux 
trapping. However, over an azimuthal angle of ~ 1°, a 
residual fringe field reduces the storage ring central field 
by about 600 ppm, increasing to 3000 ppm at the edge 
of the aperture. This fringe field reduces the field ho- 
mogeneity at large positive radial and negative vertical 
distances (Fig. |l|). 

The magnetic field is weighted with the muon distri- 
bution. Four thousand simulated muons were tracked 
for 100 turns through a field mapped by trolley measure- 
ments to evaluate the field muons encounter. The result 



agrees with B measured with the trolley probes, averaged 
over the azimuth, and taken at the beam center. 

Two largely independent analyses of ujp, which is 
proportional to (B), were made using different selec- 
tions of NMR probes. Their results agree to within 
0.03 ppm. The final value is a}p/27r = 61 791 256 ± 
25 Hz (0.4 ppm). Table || lists the systematic errors. 



TABLE I. Systematic errors for the uip analysis 



Source of errors 


Size [ppm] 


Absolute calibration of standard probe 


0.05 


Calibration of trolley probes 


0.20 


Trolley measurements of Bo 


0.10 


Interpolation with fixed probes 


0.15 


Inflector fringe field 


0.20 


Uncertainty from muon distribution 


0.12 


Others f 


0.15 


Total systematic error on ujp 


0.4 



f higher multipoles, trolley temperature and its power supply 
voltage response, and eddy currents from the kicker. 



The frequency LUa is obtained from the time distribu- 
tion of decay counts. The e"*" are detected by calorime- 
ters whose photomultiplier signals have a typical FWHM 
of 5 US. The signals are sampled every 2.5 ns by 8-bit 
waveform digitizers (WFD) with at least 16 samples per 
e+ event. The samples are used to determine 6+ times 
and energies, and for pulse overlap studies. The pulse- 
reconstruction algorithm fits signals above baseline to an 
average pulse shape, determined for each calorimeter in- 
dividually. Multiple pulses can be resolved if their sep- 
aration exceeds 3 to 5 ns. Systematic effects associated 
with the algorithm were extensively studied. 

Whereas for the 1998 data, muon decay and spin pre- 
cession (Eq. HI) was adequate to describe the e"*" time 
spectrum, the higher count rate and much larger data 
sample for 1999 required careful consideration of (1) e+ 
pulses overlapping in time (pileup), (2) coherent betatron 
oscillations, (3) beam debunching, (4) muon losses, and 
(5) detector gain stability during the muon storage time. 

(1) The number of pileup pulses in the reconstructed 
data is proportional to the instantaneous counting rate 
squared and to the minimum pulse separation time of the 
reconstruction algorithm. It amounts to 1% when the fits 
of LUa to the data are started, and distorts the e+ time 
spectrum because of misidentification of the number, en- 
ergies, and times of the positrons. Since the phase 4>a{E) 
(cf. Eq. 1^) depends on the energies of the positrons, 
pileup has a phase which differs from (j)a leading to an 
error in oja ■ Therefore, the data are corrected prior 
to fitting by subtracting a constructed pileup spectrum. 
Positrons found within a window at a fixed time after the 
pulse that triggers the WFD are treated as if they over- 
lap with the trigger pulse. The width of the window is 
taken equal to the minimum pulse separation time. Only 
data with energies at least twice the hardware threshold 
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are fully corrected with this method. Our 1 GeV hard- 
ware threshold leads to a choice oi E > 2 GeV in the oJa 
analysis. 

The contribution to pileup from signals too small to 
be reconstructed is not accounted for by the procedure 
described above. These small signals distort the pulse 
reconstruction but do not, on average, affect the en- 
ergy. However, such unseen pulses introduce small time- 
dependent shifts in (j^aiE) and A{E). The time depen- 
dence of the asymmetry, being more sensitive than the 
phase, is used to set a limit on the shift of uja ■ 

Fig. ^ shows the agreement between the positron en- 
ergy spectrum after pileup subtraction and the spectrum 
when pileup is negligible, together with the uncorrected 
spectrum. The inset illustrates that the average energy 
after pileup subtraction is constant with time, as ex- 
pected. Without accounting for pileup, a shift in uJa of 
0.3 ppm would result. 





> 

2.4- 








' — s 

>. 2.39 - 




100 200 300 

Time [|j,s] 
"^STv-^ ^^^^^ 



1 2 3 4 5 6 

E [GeV] 

FIG. 2. The energy spectrum of the detected positrons 
above 1 GeV at all times (thick line) and at only late times 
(thin line) together with the pileup-subtracted spectrum at all 
times (dashed line). The inset shows the energy above 2 GeV 
averaged over one g — 2 period as a function of time before 
(filled circles) and after (open circles) pileup subtraction for 
a typical detector. 

(2) The storage ring is a weak focusing spectrometer 
(field index n = 0.137) with an aperture which is large 
compared to the 18(w) x 57(h) mm^ inflector aperture. 
Therefore, the phase space for the betatron oscillations 
defined by the acceptance of the storage ring is not filled. 
In combination with imperfect injection angles and an 
imperfect horizontal injection kick, this results in beta- 
tron oscillations of the beam as a whole - coherent beta- 
tron oscillations (CBO). These oscillations are observed 
directly using the fiber monitors. They also modulate the 
positron time spectra, since the acceptance of a calorime- 
ter depends on the muon decay positions. The dominant 
effect is due to the horizontal oscillations, which decay 
with a time constant of ~ 100 /is. 

(3) The cyclotron period of the bunched beam leads 
to a strong modulation of the decay positron time spec- 
tra which remains at 32 /is when we begin our fits. The 
modulation structure is eliminated from the time spectra 
by uniformly randomizing the start time for each detec- 
tor and each storage fill over the range of one cyclotron 



period. 

(4) Losses of muons during the data collection are min- 
imized by controlled scraping j^] before the data collec- 
tion is started. The magnitude and time-dependence of 
remaining losses are studied using coincident signals in 
the front scintillation counters of three adjacent calorime- 
ters. 

(5) Detector gain changes and time shifts are moni- 
tored with a pulsed laser system and also using the e+ 
energy spectra. From 32 /is on, the gains of the detec- 
tors except two are stable to within 0.1% over 10 dilated 
muon lifetimes. The reconstructed times are stable to 
within 20 ps over 200 /iS (0.1 ppm). 

The raw WFD data were converted into 6+ energies 
and times using two independent implementations of the 
pulse reconstruction algorithm. Four independent anal- 
yses of Ua were made. For simplicity of presentation one 
will be described and the principal differences of the oth- 
ers will be discussed. 

The pileup corrected time spectrum for the sum of the 
detectors is fitted by the function, 

f{t)^N{t)-b{t)-l{t). (3) 

Here, N{t) is the muon decay and precession function of 
Eq. ^ and h{t) is the coherent betatron oscillation func- 
tion, 

h{t) = \ + Ai, e-*'/^' cos(wf,t + </){,), (4) 

and l{t) is the muon loss term, 

l{t) = \ + nie-''^'. (5) 

The quantities A;,, t;,, and (/>f, denote the CBO am- 
plitude, time dependence, angular frequency, and phase, 
respectively, and ni and t; denote the fraction of lost 
muons and its time dependence. The CBO frequency Uh 
is determined from a Fourier analysis of residuals in a 
fit of Eq. H to the data. The remaining 10 parameters in 
Eq. ^ are adjusted, in the sense of minimizing x^. The fre- 
quency (J correlates strongly only to (pa- Consequently, 

is insensitive, unlike fo fl^6 values of the other 8 
parameters and to the functional forms of h(t) and l{t). 
Fig. H demonstrates the good agreement of data and fit, 
as indicated by = 3818 for 3799 degrees of freedom 
(dof). ^ ^ 

The internal consistency of the results was verified in 
several ways. As an example. Fig. ^ shows the results 
when the fit range starts at increasing times after in- 
jection and UJa is seen to be constant within statistical 
errors. Fig. jija shows the results for fits to the spectra 
from individual detectors (x^/dof= 30/21). The result 
for LOa obtained from the average of individual detector 
fits (Fig. ^) is consistent with the fit to the sum (Fig. 
^) to within 0.07 ppm. The fitted lifetime, after cor- 
recting for muon losses, agrees with that expected from 
special relativity to better than 0.1%. 
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FIG. 3. Positron time spectrum overlaid with the fitted 10 
parameter function (x'^/dof= 3818/3799). The total event 
sample of 0.95 x 10^ e+ with E > 2.0 GeV is shown. 
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FIG. 4. a) The fitted frequency iJa/2TT vs. start time of fit 
is shown together with representative error bars. The band 
indicates the size of expected statistical fiuctuations. b) Fits 
to the spectra from individual detectors. Detectors 2 and 
20 were excluded from the uia analyses because of a readout 
problem and a miscalibration, respectively. 



Two of the three other analyses used fitting functions 
similar to Eq. ^. The principal differences between the 
analysis described and the others consist of a somewhat 
different choice of data selection and fitting parameters, 
a refined treatment of detector gain changes, and alterna- 
tive ways to determine pileup. In one analysis, a pileup 
correction was made by varying the minimum pulse sep- 
aration time in the reconstruction algorithm, whereas 
in the other analysis pileup is incorporated in the fit- 
ted function. The pileup phase was fixed to the value 



obtained from the constructed pileup spectrum. In the 
fourth analysis, the data are randomly split in four sam- 
ples ni - 714 which are rejoined in u{t) — ni{t) + n2{t), 
v{t) = n-i{t — Ta/1) + ni{t + Ta/1), and the ratio 

u{t) - v{t) 



r{t) 



= A{E)siXi{uJat + (j)a{E))+€ (6) 



u{t)+v{t) 

where Tq is an estimate of the g — 2 period and the con- 
stant e ^ 1. This ratio is largely independent of changes 
of observed counts on time scales larger than Tq, 6.17. the 
muon lifetime. The ratio can thus be fitted with fewer 
free parameters and its results have somewhat different 
systematic uncertainties. 

The results from the analyses are found to agree, on 
LJa to within 0.3 ppm. This is within the statistical 
variation of 0.4 ppm expected from the use of slightly 
different data in the respective analyses. We combined 
the results to uJa/2TT = 229072.8±0.3 Hz (1.3 ppm). The 
only correction applied to our result was -f 0.81±0.08 ppm 
for the effects of the electric field and vertical betatron 
oscillations 0]. The error reflects the total uncertainty, 
accounts for the strong correlations between the individ- 
ual results, and is dominated by the statistical contribu- 
tion. The systematic errors are listed in Table ^ The 
uncorrelated errors were added in quadrature, while the 
correlated errors were added linearly. Most of the system- 
atics were common to all four analyses. Spin resonances, 
fit start time, and clock synchronization were consid- 
ered and each was estimated to be less than 0.01 ppm. 

TABLE II. Systematic errors for the uja analysis. 



Source of errors 


Size [ppm] 


Pileup 


0.13 


AGS background 


0.10 


Lost muons 


0.10 


Timing shifts 


0.10 


E field and vertical betatron oscillation 


0.08 


Binning and fitting procedure 


0.07 


Coherent betatron oscillation 


0.05 


Beam debunching/randomization 


0.04 


Gain changes 


0.02 


Total systematic error on uja 


0.3 



After the u)p and LOa analyses were finalized, separately 
and independently, only then was the anomalous mag- 
netic moment evaluated. The result is 

a + = = 11 659 202(14)(6) x 10"^° (1.3 ppm) 

A — i? 

(7) 

in which R ~ LJa/u)p^ /i^ — eh{\ + a^)/(2m^c), and 
A = ^i^/^ip = 3.183 345 39(10) |]. This new result is in 
good agreement with previous measurements 0,^,^ and 
reduces the combined error by a factor of about three. 

The theoretical value of in the standard model (SM) 
pO[ has its dominant contribution from quantum elec- 
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trodynamics but the weak and strong interactions con- /r^Trr^^^N , ... n / 100 GeV\^ ^ 

tribute as well. The value "''(^USY) « 140 x IQ-" (^^— j tan/3. (10) 

a^(SM) = 11 659 159.6(6.7) x lO^i" (0.57 ppm) (8) jf ^g^^i^e the difference a^(exp) - a^(SM) to 



is the sum of a^(QED) = 11 658 470.56(0.29) x 10^1° 
(0.025 ppm), a^(weak) = 15.1(0.4) x lO^^" (0.03 ppm), 
and a^(had) = 673.9(6.7) x 10-'^° (0.57 ppm). The 
term a^(QED) is obtained using the value of a from 
ae(exp) = ae(SM) |10|, and terms through order are 
included. The term a^(weak) includes electroweak con- 
tributions of up to two- loop order. The term a^(had) 
arises from virtual hadronic contributions to the photon 
propagator in 4*'* order and 6*'' order, where the latter 
includes light-by-light scattering. The 4*'* order term, 
which provides the largest contribution and uncertainty 
to a^(had) is obtained from measured hadron production 
cross sections in e+e~ collisions and hadronic r decay 
^t| . Additional data on e+e" collisions from Novosi- 
birsk a nd from Beijing ||l^ and on r decay from 
Cornell||l4[| have not yet been included in the evaluation 
of a^(had). 

In Fig. 0, the five most recent measurements of are 
shown along with the standard model prediction. The 
difference between the weighted mean of the experimen- 
tal results, a;,(exp) = 11 659 203(15) x lO'^" (1.3 ppm), 
and the theoretical value from the standard model is 



,(exp) - a^(SM) = 43(16) x 10 



-10 



(9) 



The error is the addition in quadrature of experimental 
and theoretical uncertainties. The difference is 2.6 times 
the stated error. 
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FIG. 5. Measurements of and the standard model pre- 
diction with their total errors. 

Many speculative theories predict deviations from the 
standard model value for These include super- 

symmetry, muon substructure, and anomalous W cou- 
plings. The muon anomalous g value is particularly 
sensitive to supersymmetry whose contributions to 
a^i come from smuon-neutralino and sneutrino-chargino 
loops. In the limit of large tan/3, which is the ratio of the 
vacuum expectation values of two Higgs doublets, and for 
a degenerate spectrum of superparticles with mass m, 



a^(SUSY), for tan/5 in the range 4 - 40, then m w 120 
- 400 GeV. 

In 2000, approximately four times the total number of 
positrons were recorded as compared to our 1999 data. 
Measurements are now underway with negative muons, 
which will provide a sensitive test of CPT violation. 
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